Global well-posedness for a Smoluchowski equation coupled 
with Navier-Stokes equations in 2D. 

p. Constantin 

Department of Mathematics, The University of Chicago 
5734 S. University Avenue, Chicago, II 60637 
email const@math.uchicago.edu 
and 

Nader Masmoudi 
Courant Institute, New York University 
251 Mercer St, New York NY 10012 
email:masmoudi@cims. nyu.edu 

February 2, 2008 



Abstract 

We prove global existence for a nonlinear Smoluchowski equation (a nonlinear Fokker- 
Planck equation ) coupled with Navier-Stokes equations in 2d. The proof uses a deteriorating 
regularity estimate in the spirit of (see also fl| ) 
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1 Introduction 

Systems coupling fluids and particles are of great interest in many branches of applied physics 
and chemistry. The equations attempt to describe the behavior of complex mixtures of parti- 
cles and fluids, and as such, they present numerous challenges, simultaneously at three levels: 
at the level of their derivation, the level of their numerical simulation and that of their math- 
ematical treatment. In this paper we concentrate solely on one aspect of the mathematical 
treatment, the regularity of solutions. The particles in the system are described by a probabil- 
ity distribution f{t,x,m) that depends on time t, macroscopic variable x G M", and particle 
configuration m € M. Here M is a smooth compact Riemannian manifold without boundary. 
The particles are transported by a fluid, agitated by thermal noise, and interact among them- 
selves. This is reflected in a kinetic equation for the evolution of the probability distribution 
of the particles ([21 E])- The interaction between particles - a micro-micro interaction - is 
modeled in a mean-field fashion by a potential that represents the tendency of particles to 
favor certain coherent configurations. The interaction between particles occurs only when the 
concentration of particles is sufficiently high. Mathematically, this term is responsible for the 
nonlinearity of the Smoluchowski (Fokker-Planck) equation, and physically, it is responsible 
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for nematic phase transitions. Because the particles are considerably small, and for smooth 
flows, the Lagrangian transport of the particles is modeled using a Taylor expansion of the 
velocity field. This gives rise to a drift term in the Smoluchowski equation that depends on the 
spatial gradient of velocity. It is a macro-micro term, and it causes mathematical difficulties 
in the regularity theory. 

The fluid is described by the incompressible Navier-Stokes equations. The microscopic 
particles add stresses to the fluid. This is the micro-macro interaction and it is the most 
puzzling and important physical aspect of the problem. Indeed, while a macro-micro inter- 
action can be derived, in principle, by assuming that the macroscopic entities vary little on 
the scale of the microscopic ones, the "scaling up" of the effect of microscopic quantities to 
the macroscopic level is more mysterious. A principle based on an energy dissipation balance, 
and that recovers familiar results in simple cases was proposed in [6], where the regularity 
of nonlinear Fokker-Planck systems coupled with Stokes equations in 3D was also proved. 
The linear Fokker-Planck system coupled with Stokes equations was considered in [19]. The 
nonlinear Fokker-Planck equation driven by a time averaged Navier-Stokes system in 2D was 
studied in [7j. 

An approximate closure of the linear Fokker-Planck equation reduces the description to 
closed viscoelastic equations for the added stresses themselves. This leads to well-known non- 
Newtonian fluid models that have been studied extensively. For regularity results we refer 
to Lions and Masmoudi [18] where the existence of global weak solutions was proved for an 
Oldroyd-type model. In Guillope and Saut [13] and [H], the existence of local strong solution 
was proved. Also, Fernandez-Cara, Guillen and Ortega [TT], [TO] and [12] proved local well 
posedness in Sobolev spaces. We also mention Lin, Liu and Zhang [16] where a formulation 
based on the deformation tensor is used to study the Oldroyd-B model. 

An other model for the polymers is the FENE dumbbell model. From mathematical point 
of view, this model was studied by several authors. In particular W. E, Li and Zhang [9], 
Jourdain, Lelievre and Le Bris |jl5j and Zhang and Zhang ,20j proved local well-posedness. 
Moreover, Lin, Liu and Zhang [T7] proved global existence near equilibrium. 

1.1 The model 

Consider the system 

^ + vVv- vAv + Vp = V • r in S7 x (0, T) 

< |£ + ^;.v/ + div,(G(t;,/)/)-Ag/ = in l^x(0,T) (1) 

di-vv = in Ox(0,r), 

where Tij = Jj^j 7lj\rn)f{t, x, m)dm + f^j Jj^j ^^^\mi,m2)f{t, x, mi)f{t, x, m2)dm. We denote 
G{v, f) = V gU + W where W = cidjVi and U = Kf is a potential given by 

Uit,x,m)= [ Kim,q)f{t,x,q)dq (2) 

with a kernel K which is a smooth, time and space independent symmetric function K : 
M X M ^ R. We also take = M^. 
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1.2 Statement of the result 



Theorem 1.1 Take v{0) € M^i+^0''"nL2(M2) and /(O) G W^'^'^R-''), for some r > 2 and Eq > 
and f> 0, Jj^ /o G n L~. Then (Cp has a global solution in v e L^^{W^''') n L^^iW^''') 
and f £ L~ Moreover, for T > Tq > 0, we have v G L~((ro,r); VF^-^'O- 

1.3 Prehminaries 

We define C to be the ring of center 0, of small radius 1/2 and great radius 2. There exist 
two nonnegative radial functions x ^^i^d f belonging respectively to 'D{B{0, 1)) and to T>{C) 
so that 

q>0 

\p-q\>2^ Supp V9(2-^-) n Supp V3(2-P-) = 0. (4) 
For instance, one can take x ^ ^(-^(0, 1)) such that x = 1 on -6(0, 1/2) and take 

9^(0 = x(e/2)-x(e). 

Then, we are able to define the Littlewood-Paley decomposition. Let us denote by the 
Fourier transform on R'^. Let h, h, Aq, Sq {q G Z) be defined as follows: 



h = T ^Lp and h = ^x-, 
Aqu = T-Hip{2-'^0^u) = 2^' J h{2'^y)u{x - y)dy, 

Squ = T-\x{2-'0^u) = 2'^^ j h{2'iy)u{x - y)dy. 

We use the para-product decomposition of Bony ([3]) 

uv = T^v + TyU + R{u, v) 

where 

TuV = ''^^Sq-luAqV and R{u,v) = Aq'uAqV. 

qei k-g'|<i 
We define the inhomogeneous and homogeneous Besov spaces by 

Definition 1.2 Let s he a real number, p and r two real numbers greater than 1. Then we 
define the following norm 



def, 



W\b= - ||'S'o'"||lp + 



i2'^'\\AqU\\Lp) 



and the following semi-norm 



def 



i2^^\\AqUU.)q^^ 



r(z) 



Definition 1.3 



• Let s be a real number, p and r two real numbers greater than 1 . We denote by ^ the 
space of tempered distributions u such that \\u\\^^ is finite. 
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If s < d/p or s = d/p and r = 1 we define the homogeneous Besov space Bp ,, as the 
closure of compactly supported smooth functions for the norm \\ ■ 



p,r 



We refer to [4] for the proof of the following results and for the multiplication law in Besov 
spaces. 

Lemma 1.4 

W^quhb < 2'^^a-T;)i\\Agu\\L'' for 6 > a > 1 
The following corollary is straightforward. 

Corollary 1.5 Ifb>a>l, then, we have the following continuous embeddings 

S \ab 

Definition 1.6 Let p be in [1, oo] and r in M; the space Ll^{C^) is the space of distributions u 
such that 

def 

\Mlp(ot-c-) - sup2«''||Aynj|iP(ioc) < oo. 

' ' ' q 

We will use the following theorem from [5] 

Theorem 1.7 Let v be the solution in Lj,{H^) of the two dimensional Navier-Stokes system 

dv 

i^S-X dWv = 

v\t=0 = Vq. 



„ + V ■ Vv — uAv = —Vp + / 
ot 



with an initial data in and an external force f in L}p{C ^) R L'^{H then, for any e, 
a Tq in the interval ]0,T[ exists such that 

2 A deteriorating regularity estimate 

The main part of this section is the proof of a deteriorating regularity estimate for transport 
equations in the spirit of [1] and [5]. After this proof, we will apply this estimate in order to 
prove Theorem I l.li 

We also denote H = {-Ag + 1^/^ with s> d/2 + 1. 

Theorem 2.1 Let a and (5 he two elements o/]0, 1[ such that a + (3 <\. A constant C exists 
that satisfies the following properties. Let T and A be two positive numbers and v a smooth 
divergence free vector field so that 

a — v\\-^-i^f^^o-^ > (3 and a + A|| VfUji^j-^o^ < 1 — /3. (5) 
Consider two smooth functions f and v so that f is the solution of 

(6) 



dtf + v-Vf + divgiG{v,f)f)-Agf = 



f\t=o — /( 



0- 
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Then we have, if X > 3C, 

MIU) < 3||/obj^^^(^,-.) + —Aq+\v) (7) 



where 



Ml{v) =^ sup 29'^-*'J.^W||Agi>(t)||iP or (8) 

iG[0,T],g 

Ml{f) =^ sup 2^'^-*-^W||A,/(t)||i,(^-.) w^th (9) 

iG[0,T],g 



Jt' 



(10) 



def 

We will use the notation fq = Ag/. Applying the operator Ag to the transport equation 
dS]), we get 

dtU + Vi^ • V/, + div^(G(5g_iT;, 5,-i/)/g) - Ag/g + iig(t;, /) = ^^^^ 



/?|t=o - N/o- 

•est term. 

We denote 



where Rq is a rest term. 



Nl(t,x) = [ \Hfq\^dm (12) 
Applying H to (jlip and taking the norm on M, we get 

dtN^ + Sq.iv ■ VNl + V{Sq.iv, A) + iVg^/gl' + / Hfq{HRq{v, f))dm = (13) 

J A/ 



where 



V{v,h,f)=d,Vi [ {Hd\vg{4f)){Hf)dm.+ f {Hd\Yg{Vghf)){Hf)dm. (14) 

J A/ JM 

Hence, arguing as in [7], we have |y(5g_if , /g)| < C(|V5g_if| + | l^g-i/l |l2(m))A''2. 
We will use now the following lemma, postponing its proof: 

Lemma 2.2 Rq{v,f) satisfies 

2'''^~'^.Mt)\\HRqiv{t),fm\LnL-) < Ce""'"^^"^^'^") 

'M^+\v) + ^l + \\SqVv{t)\\L^+ l|Ag.Vt;(t)|U^)M,-(/)y (15) 

Taking the LP norm of Nq^ we get 



||iV,(t)||Lp < ||iVg(0)||LP + / ||//i?g(7;(t'), /(*')) + ||VV(i')llL-||iV,(t')llLpdt'. 

JO 

After multiplication by l^^^^i-^-^^) ^ we get 
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Jo 
Jo 

Then, using the inequahty p3]) and taking the sup over q, we get 

(16) 

M^+Hv) + M^if)(^l + 2\\SgVv{t')\\L^+ \\^,''^vit')\\L^)]dt'. (17) 

\q'-q\<N ^ 

As AllVt^lIji ((^0) is smaller than (a — /3), we have 

^CA||V«||ji^(pO) ^ gC(<7-/3)_ 

Moreover, by definition of ^q^\{t^t'), it is obvious that 



te[0,T],gJo 



/■*2-*,.A{t/)(||5 Vz;(tOl|Loo + \)di! < -i 



Then, we obtain that 



'g2 



MK/) < ||/obj^,^(j/-) + i^M^"+^(«) + C||VHIz^(co)Mr(/) + i^Ma (/) 
This proves the theorem of course if we prove the estimate (|15|) of the lemma. First of all, let 
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us decompose the operator Rq. We have 

6 

Rii^J) = Y^RlivJ) with 

i=i 

d 
d 
d 
d 

R>^f) = i;div,(c^JA,(T^a,V)) + div,(A,(T^V5C/)), 

^.'(^'/) = i;div,(4^[A„r9^.,.]/) + div,([A„rv,d/) 

+ di^s (^(VpC^, /) + A,_i VgC/A,+i/, - A,_2Vgi7A,_i/,) 

Indeed, 

d 

2 d 

= ^ /) + 5; T,,a, A,/ + Agi?!^;^ 5,/), 
^=1 j=i 

Then, we use that 

d 

j=l \q-q'\<i 

Iq-9'!<i 

= Sq-iV^djfq + Aq_iV^djAq+ifq - Aq_2V^ djAq_ifq 

Hence, 

3 

• V/) = 5] R'qiv, f) + Sq^.V ■ Vfq. 

In the same way, we have 

6 

A,(div,(G(^;, /)/)) = J2 /) + divg{G{Sq.,v, Vi/)/?)- 

e=4 



Let us estimate the six terms appearing above. We have 

Let us begin with Rg{v, /). By definition of the paraproduct, we have 



d 

j=l q' 



As, if |g — gr'l > 2 then the above term is equal to 0, we deduce that 

\\HRl{v{t)J{t))h,^L2^<C Yl ll^5,,_iV/|Uoc(i2)||A,,^(t)|Up. 

\q-q'[<2 

Using the fact that, if \q - q'\ < 2, then \\HSg>^iV f\\Loo^L2) < C2«||iJ/(t) 11^00(^2) < C2i, we 
infer that 

\\HRl{v{t)jm\LriL^)<C2'^ E \\^q'vmL^<C E \\VAg,v{t)h.. 

k-9'l<2 \q-q'\<2 

Hence 

2?^-*.,A(*)||i7i?i(^;(t),/(t)) 11^,(^2) < CM^+\v) E 2-^fo\\s.^^in\\Loodt'+xf*\\s^,Vv{t')Uoodt'_ 

k-9'l<2 

But, it is obvious that 

/* \\Sg,Vv{t')\\Loodt' - f \\SqVv{t')\\Loodt' < f 11(5,, - Sg)Vv{t')\\Loodt'. 

Jo Jo Jo 

Using the fact that \q — q'\ < 2, we get 

\\Sg,Vv{t')\\L00dt' - \\SgVv{t')\\L00dt' < C \q " | 1 1 1 1 (^Q ) • (18) 

So it turns out that 

2'^-^^Mt)\lHRlivit),fm\LnL^) < 2'''^^'''Hio°)M-+\v). (19) 
Now let us look at Rq{v, /). By definition of the paraproduct, we have 

d 

j=l q' 

d 

j=l q' 

The terms of the above sum are equal to except if |g — g'| < 2. Moreover, by definition of 
the operators A,, we have 

[Sg,_iv^,Ag]djAg,f{x) = 2'^^ [ /i(2^(x " v)) {S g> _iv^ (x) " S g^ _iv^ (v)) dj Ag, f {y)dy . 

jRd 

So we infer that 

\\H[Sg>_,V^,Ag]^,Ag,f{x)h2^M)<2-^VSg>_,V\\L^2''^(^(^^^^^^ 



Hence, 

\\H[Sg>.iv^,Ag]d,Ag>f{x)h.iLHM))<2-'\\VS,>^iv\\L^\\Hd,^ 
Then, we have, using Inequahty ([TS]) . 

<CM^if) 2''"''"2^(^^Hl|V(5,'-i-5,)t;(t)|U^ + ||V(*)llL-)- 

k-'?'l<2 

So, we get 

2<i'^-'^.At)\lHRlivit),fit))hr.^L^) < CM,-(/)2^""''l'£^(o^) 

)\SgVv{t)\\L^ + ^ ||V(A,,7;(t)||Loo). (20) 



lg-9'l<2 



For we have 



\q'-q"\<l 
q'>q-2 



< C 2'^-'''\\A,,Vv\\Lp\\Hfh. 

q'>q-2 



Hence, 



q'>q-2 

Then, we see that the sum converges since 

l^,,A(t) - '^q'At)\ < A||vH|j^(^o)k - <i'\ <{<y-m- q'\ (21) 

and l + cr — (cT — /?) = ! + /?> 0. Hence, we get 

2'i^-'^^--^)\\HRl{v,f)U,(L.^ < CMl-^\v)\\Hf\\L^(L.y 

The estimate for Rq{v, f) = Rq'^{v, /) + Rq''^{f) is the same as the estimate for Rq{v, /). 
Indeed, we have 

\\HR'/{v{t),fm\L,^L^) < C Y \\'^9HSq,^if\\Loo^L2)\\Aq,Vv{t)\\Lr 

k-9'l<2 

< C Y l|A,'V^;(t)||LP 
k-'?'l<2 

where we used that || Vgi?S'g'_i/||/^oo(^2) < C. Hence, we conclude as for Rq{v,f). Besides, 

k-g'l<2 

< c Y iW'fimL. 

\q-q'\<2 



Hence, we conclude as for Rl{v, f) and get 

2^--*-^W||ifi?f < 2^'^^^''H(o'>^M^{f). (22) 

We write Rg{v, f) = Rq'^{v, f) + Rq''^{f). The estimate for Rq{v, /) is similar to the one 
for Rq(v, f) with the only difference that we have to use the regularity of Vv. We have 

d 

i = l q' 

The terms of the above sum are equal to except if \q — q'\ < 2. Moreover, by definition of 
the operators Ag, we have 

[Sq:-idjv\ Aq]Aq>f{x) = 2'^^ [ /i(2'?(x - y)){Sq: .id^ (x) - Sq> ^id,v\y)) Aq, f {y)dy . 

So we infer that 

\\HRl'\v,f)\\L2^M) < 2-'^|v25,._iH2^'^((2^| • I X |/i(2^.)l) * || V,i?A,,/|U2(M)) (re). 
Hence, 

||lfi?^.l(z;,/)|j^,(^2(M)) <2-''||V25,,_it;||ip||V,i/A,,/||i.o(i2(M)). 

Then, we have, using Inequality ([H]), 

\q~q'\<2 
q"<q'~l 

Hence, 

q"<q+l 

and the sum is uniformly bounded since a — 1 + X\\Vv\\j^i ^^o) ^ —P- Then, we argue in a 
similar way for H-f/'-Rg ^(/) ||lp(l2{m)) and get 

and we conclude as above with M^^^{v) replaced by M^{f). 

Finally, the estimate for Rq{v, /) is exactly the same as the one for Rq{v, f) since, we also 
have that || Vg-?//||^cx)(^2) < C. 

3 Global existence 

Now, we turn to the proof of our main theorem. First, we notice that the local existence in 
with V G L;~([0,r);T^i''^) nL2,^([0,r);1^2,r) j ^ Lf^^{[{),T);W^^'{H~')) can be easily 

deduced from standard arguments. Moreover, from regularity estimates for the heat equation, 
we have for all < Tq < T, v € ^^^((To, T); VF^-e.r)^ 
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We want to prove that we can extend the solution beyond the time T. It is enough to 
prove that S/v G L°°((0,r) x M^). 

The local existence result tells that, for any Tq in ]0, T[, the solution {v,f) of ^ belongs 
to the space L^^([To, r[; T4^2-£,r ^ w^^''{H-')) for any e > 0. Sobolev type embeddings of 
Corollary [T3] imply that 

Choosing s < 1 — 2/r and p = oo in the above assertion implies that {v,t) E L'^^{C^~^'^ x 
C"{H~'^)) where a = \ — e — 2/r > 0. So we can apply Theorem 11.71 and we can choose Tq 
such that, with the notations of Theorem 12.11 we have 



min{a — P,l — a — P) 

" "^[To,T,(^") - 3A 

The deteriorating regularity estimate of Theorem 12.11 applied with a and between Tq and T 
tells exactly that / satisfies 

MUf) < 3||/||c^(H-.) + —M^+Hv). (23) 

Now, we have to estimate Vf . The two dimensional Navier-Stokes equation can be written as 

dtv - I'Av = P{v ■ Vv) + PDt 

where P denotes the Leray projector on the divergence free vector field. Exactly along the 
same lines as in the proof of Theorer d2.1l we have 

2«("+i)-*'.^W||P(fVf)-P(5qZ;-VAgt;)||Lco < CM^+\v)(\\Sg\7v{t)\\Lo^ + Y,'^'"''\\'^\'vit)\\L 

^ q'>q 

Moreover, it is obvious that 

2.(<x-i)-*,..(t)||p(5^^ . vA,^)||ioo < C\\v{t)\\^iM^+\v). 
So it turns out that 

29(-+i)-'f,.AW||A^P(^;.Vz;)||ioo < CM^+\v)(\\SgVvit)\\Lo^ + Y, 2^''-'''^\\Vv{t)\\L^+2'^\\v{t)\\^ 

^ q'>q 

(24) 

Using well known estimates on the heat equation (see for instance ^) and Inequalities (I23p 
and dMl) , we get that 

M^+\v) < \\vo\\c^+i + (^j + 2'^Fg{To,T)^Mliv) + ^M^ir) 



with 



F,{n,T)'M sup re-2^^M||,(t')|l .dt'. 
t&[To,T]JTo 



Holder inequality implies immediately that 



F,(To,r)<-^2-^||^;|| . 

z/4 ^Tf,.T]\" ) 
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Moreover, it is easy to see that 



So, we infer that 



3C„ „ fC C C 

V 



M^+\v) < \\vo\\c^+i + — ll-^ollc- + ( T + — + —11^11.4 



Now it is enough to choose To such that the quantity 

'C C _C 

3 



is small enough. Then as a is greater than 0, the solution {v, r) of the system ([T]) is such 
that (Vf ,r) belongs to L°°([To,T] x M^); this concludes the proof of Theorem I l.li 
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